Duality and fluctuation relations for statistics of currents on cyclic graphs 
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We consider stochastic motion of a particle on a cyclic graph with arbitrarily periodic time dependent kinetic 
rates. We demonstrate duality relations for statistics of currents in this model and in its continuous version of a 
diffusion in one dimension. Our duality relations are valid beyond detailed balance constraints and lead to exact 
expressions that relate statistics of currents induced by dual driving protocols. We also show that previously 
known no-pumping theorems and some of the fluctuation relations, when they are applied to cyclic graphs or to 
one dimensional diffusion, are special consequences of our duality. 
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I. INTRODUCTION 



Dynamics of mesoscopic structures can often be modeled 
as a stochastic motion on a network (H-0]. Nodes of this 
network represent metastable states of a structure and links 
correspond to allowed transitions. Examples can be found 
among biological enzymes JHI^l], electronic circuits 10,111] and 
heat pumps [9]. Typically, nanomechanical systems should be 
able to perform cyclic operations, which correspond to motion 
along the cycles in such networks JUS]. 



In this article, our main focus will be on stochastic particle 
motion in a cyclic graph. In such a graph, nodes are connected 
one by one, forming a closed chain, i.e. a cycle. Number of 
nodes of a cyclic graph is equal to the number of its links and 
each node has exactly two neighboring nodes to which it is 
connected. This graph is already a sufficiently general repre- 
sentation of the kinetics for many nanoscale structures. Un- 
derstanding dynamics on such a graph is also a starting point 
for understanding more complicated systems. For example, 
arbitrary networks can be decomposed into fundamental cy- 
cles In such a graph, we will study counting statistics of 
current induced by time-dependent driving and by forces that 
explicitly break the time-reversal symmetry. 



Our main result is that the Markovian kinetic on cyclic 
graphs with arbitrary time-dependent kinetic rates has exact 
duality for statistics of particle currents. This duality leads to 
new exact relations for statistics of currents in systems driven 
by time-dependent protocols. These relations are akin to 
known fluctuation theorems for currents [10] but they are ap- 
plied to systems with arbitrary time dependence of all kinetic 
rates. In fact, for cyclic graphs, we show that known fluc- 
tuation relations, as well as other known exact results called 
no-pumping theorems, are special consequences of our dual- 
ity relations. 



H. MASTER EQUATION AND ITS DUAL EQUATION ON A 
CYCLIC GRAPH 

Consider a cyclic graph with N nodes and N links. Links 
represent allowed transitions between nodes. We assume that 
all links can be traversed in both directions with some rates. 
Detailed balance is not assumed. Let kf and fcr be the kinetic 
rates of transitions from state i into, respectively, states i + 1 
and i — 1, as shown in Fig. 02a). We will identify indices i = 
as i = N, and i = N + l as 2 = 1 if necessarily. The evolution 
of the probability vector p is given by the master equation: 



d t p = Hp, 



(1) 



where we will call H the master operator. Let (j\ and \i) be 
the bra- and ket- iV-vectors with the only nonzero unit compo- 
nents at j-th and i-th positions respectively. In this basis, the 
operator H has components Hi±\^ = kf, H. L ,i = —k~^ — k^, 
and zero otherwise. The constraint of column-sum-zero on 
the master operator, 



N 



(2) 



i=i 



implies the conservation of the probability on a graph, 

L,i=iPi = 1 - 

The master operator, H, can be written as a product of two 
operators: 



H=QJ, 



(3) 



where 



Q\i) = \i - 1) - \i) J\i) = kr\i)-kf\i + l). (4) 

We will need the following property, namely, with each master 
operator H on a cyclic graph, we can associate a dual master 
operator, given by full : 



Q T J T , 



(5) 
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where T is a transposition operation. In components, we have 

The operator, H d , satisfies the probability conservation condi- 
tion Eq. d2j and hence it is a valid master operator describing 
stochastic particle motion on a cyclic graph with kinetic rates: 

kf + = K, kf- = k+_ v (7) 

m. STATISTICS OF CURRENTS ON A CYCLIC GRAPH 

Suppose that J is a random variable, which we call the 
current, defined so that whenever the particle makes a jump 
(transition) in the clockwise direction we increase J by unity 
(J — > J + 1), and whenever the particle makes a jump 
in the counter-clockwise direction we decrease J by unity 
(J — > J — 1). Here we assume that all kinetic rates are 
time-dependent and change according to some periodic pro- 
tocol with period t. Such a steadily driven system eventually 
enters a regime with periodically changing population prob- 
ability vector, i.e. pit + r) = pit). We also introduce the 
current per period of the driving protocol in the limit of large 
observation time, which is also called the conserved current, 

1 f nT 

J c = lim - / dtJit), (8) 

n->oo n J Q 

where n is the number of cycles of the driving protocol. In 
general, the distribution of conserved current has an asymp- 
totic form P(J dtJ(t)) ~ exp(n5(J c )) in the limit of large 
n, or written as 

P{J C ) ~ e s{r \ (9) 

where S (J c ) is called the large deviation function (LDF) lfl2ll . 
The cumulant generating function (CGF), p, x , generates cu- 
mulants of the conserved current and is defined through the 
relation (exp(x J ™ T dtJ)) ~ exp(n[i x ), or say: 

(e xJC )-e^. (10) 

The ith order derivative of /z x with respect to the auxiliary 
counting parameter x, at x = gives the ith order cumu- 
lant of the conversed current. Moreover, the CGF is con- 
nected with the LDF through the Legendre transform fj, x = 
Maxjc [x J c + S{J C )] d]. Thus, once we find the CGF, we 
are able to obtain the LDF through the inverse Legendre trans- 
form so that the distribution of the conversed current can be 
reconstructed. 

There is a well developed and reviewed approach to calcu- 
late n x for models of stochastic evolution Il4 l8l llOl 1 1 44 - 4 1 611 . 
Let Pi,j(t) be the probability that at time t the particle will be 
on the i-th node and by this time the produced current will be 
J. Then the master equation for probabilities Pi,j{t) reads 

-r;Pi,J = -(K + K) p i,J + K+ iP+iJ+i + K- iPi-i,J-x- 



Let Z = Z(x, t) be the vector of generating functions with 
components Z t = XX-oo Pi,j(t)e xJ - Multiplying Eq. CCQ) 
by e xJ and summing over J, we find that the generating func- 
tion Z satisfies the following equation: 

dtZ = H X Z, (12) 

which we will call the twisted master equation. We will call 
H x the twisted master operator. It has off-diagonal com- 
ponents given by \H X ) = Hi±i.ie ±x . The term 

"twisted" is justified here because the operator, H x , can be 
obtained from the master operator, H, in Eq. (Q]i by multiply- 
ing (twisting) its off-diagonal matrix elements, responsible for 
clockwise or counter-clockwise transitions, with auxiliary pa- 
rameter e x or e~ x respectively. 

A formal solution of Eq. (flZb can be written in the form 

Z( X ,r) = e^=U x p(Q), U x = fe£ A *®, (13) 

where U x — U x {t) is the corresponding evolution operator 
and p(0) is the vector of the initial population probabilities of 
the nodes. To obtain statistics of the conserved current, J c , we 
note that in the limit of large observation time the information 
about the initial state vector, p(0), becomes irrelevant and the 
evolution of Z is dominated by the largest eigenvalue of the 
evolution operator, U x . The logarithm of this eigenvalue gives 
the CGF of currents [i x , namely, 

U x \u x )=e^\u x ), (14) 

where \u x ) is the eigenstate of U x that corresponds to the 
largest eigenvalue. 

Similarly, for the dual master operator, H d , we can build 
its twisted version, H d , by multiplying its off-diagonal ma- 
trix elements, which are responsible for clockwise or counter- 
clockwise transitions, with e x or e~ x , respectively. The anal- 
ogous evolution operator and the CGF for the dual graph are 
given by 

U$=feK 6 & t \ U«\u x ) = e^\u x ), (15) 

with \u x ) the eigenstate of U x corresponding to the largest 
eigenvalue. 

IV. DUALITY IN TWISTED MASTER EQUATION 

The duality on the level of master operators has been dis- 
cussed previously |[Tll[T7ll . Here we will show that a similar 
duality can be introduced on the level of twisted master oper- 
ators, which describe statistics of currents in original and dual 
models. 

Let us introduce two twisted operators in cyclic graphs with 
arbitrary driving, namely: 

2x1*) = I* - !) - e x N), J x \i) = e~*kr\i) - k+\i + 1), 

(16) 
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FIG. 1: Sketch maps for parameters of the twisted master operators 
on links connected to state i for (a) original kinetic model and (b) its 
dual graph. 



FIG. 2: The duality and hidden symmetry relations among the 
four twisted discrete operators as well as their continuous correspon- 
dences, the four twisted Fokker-Planck operators. 



where the former is a constant operator while the latter con- 
tains the time-dependent transition rates. With these opera- 
tors, two twisted master operators can be built, namely 

H X = Q X J X , and i? d x = QTjJ. (17) 

One can check that H x is a twisted master operator that corre- 
sponds to the evolution of the counting statistics in the original 
kinetic model, as seen in Eq. dT2l . While if_ x is the twisted 
master operator that corresponds to the master operator H A 
of the dual model, however, with the counting parameter \ 
(— x) replaced by — \ (x)- The master operators in the origi- 
nal and the dual models, H x and H_ x , are further connected 
with each other by the duality relation for their kinetic rates: 

kf+(t) = ki(t), fcf~(t) = fctiW, (18) 

as illustrated in Fig.Q] 

Let us denote H- x = H x — J X Q X - Since H- x and 
H x differ merely by a transposition, their evolution operators 
U- x = Texp^dtH^ and U x = Texp^dtH^ 
have a "time-reversal" correspondence as a consequence of 
the time ordering operator T, namely, 

U x ,f = (U- X , B ) T ■ (19) 

Indexes "F" and "B" correspond to "forward" and "back- 
ward" protocols, such that the time-dependent parameters are 
related by k i B (t) = k iF (r — t). (In what follows, we omit 
this index if the quantities at both sides of an equation be- 
long to the same "forward" or "backward" protocols.) Sim- 
ilarly, let = {Ht x ) T = J X Q X be the "time-reversal" 
counterpart of the dual twisted master operator, H_ x - Corre- 
sponding evolution operators, U x = Texp ^ JJ" dtH x ^j and 

U d x = f exp ( J Q T dtHt^j , satisfy the relation: 

^,f=(^ x ,b) T - (20) 



H x and H x satisfy the obvious relations H X Q X = Q X H X 
and H d _ x Q x ? = Q X H_ X , and can be regarded as the super- 
symmetric counterparts to each other rfTTl I2TI1 , as well as H^ x 

and H- x . The relations among the four twisted operators are 
summarized in Fig. [2] 

These four twisted discrete operators are isospectral for 
arbitrary time-dependent protocols. This follows from the 
fact that the operator Q x is time-independent, and if \u x (t)} 
satisfies the equation d t \u x (t)) — H x \u x (t)) then the vec- 
tor \u x (t)) = Q x \u x (t)} immediately satisfies the original 
twisted master equation d t \ u x (t)} = H x \u x (t)). Since solu- 
tions of evolution equations with H x (t) and H x (t) coincide 
up to multiplication by a constant matrix, the eigenvalues of 
the corresponding evolution operators U x (t) and U x (t), are 
the same. Further, by noticing that the transposition only 
reverses the time ordering, while keeps the eigenvalues un- 
changed, the matrices U x .f, U_ x B' ^x-F an d U_ x B share 
the identical eigen-spectrum. 

As a consequence, evolution operators for counting statis- 
tics in the original model and its dual counterpart with a time- 
reversed driving protocol have the same largest eigenvalues 
and, hence, the same CGFs. We are now in a position to for- 
mulate our findings in the form of the duality relation for cur- 
rents, stating that: 

For any driving protocol on a cyclic graph with periodically 
time-dependent rates, k i (t), there is a dual protocol, defined 
as 

kf^(t) = KMr-t), kfp(t) = k+_ l>F (T-t), (21) 

with i = 1,...N. The cumulant generating functions for the 
original and the dual protocols satisfy the relation: 

Mx.f = V- X ,B- ( 22 ) 

After the inverse Legendre transformation, Eq. \22\ leads to 
the following relation between LD-functions and probability 
distributions for conserved currents generated by original and 
dual protocols, 

S F (J C )=S*(-J C ), P F (J C ) = P£(-J c ). (23) 
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Eq. d23l ) can be called a fluctuation relation in a sense that 
it describes the full statistics of currents in arbitrarily driven 
system. Eq. ( 1231 is the main result of this article. It shows 
that exact fluctuation relations can follow from duality rather 
than from relations among work, free energy and entropy. 



DUALITY AND PREVIOUSLY KNOWN EXACT 
RESULTS 



Comparing with Eq. (1211 . we have the duality relations for the 
free energy inputs: 

/&(*) = /<7f(^ - *)> A?b(*) = ft-i?(? - *)■ ( 26 ) 

The increased entropy after one clockwise cycle in the dual 
graph is exactly the decreased entropy after one clockwise cy- 
cle in the original one: 

A = J2 Hkt/K) = E Hkf + /kt). (27) 



Beyond the fluctuation relation (1231 . the duality provides 
us with a method to generate new fluctuation relations from 
already known ones. Here we demonstrate how some of the 
nontrivial known results could be derived using our duality 
transformation. 



A. No-Pumping Theorems 

When the discrete kinetic model corresponds to a system 
with strong coupling to thermal environment at inverse tem- 
perature f3, its kinetic rates, following the Arrhenius law, 
may be parameterized as kf_ x — ke^ Ei - x ~ Wi > and fcr = 
kgPlEi-Wi) w j tn some parameters k, Wi and Ei, where Ei 
corresponds to the size of the energy well of node i, and Wi 
corresponds to the size of the barrier separating nodes i and 
i — 1. In this parametrization, recalling the relation Eq. (EH . 
original and dual protocols are related as follows, 

^b(*) = -WiA-r ~ t), W t d +1<B (t) = -E itF (r - t), 

(24) 

i.e., the dual model is different from the original one by ex- 
change of energy wells {Ei} with inverted barriers {Wi}. The 
duality relation means, in particular, that if we prove some 
general property for a model in which only barriers are driven 
at arbitrary values then automatically the same property fol- 
lows for any model with the same geometry in which barriers 
are fixed but energies are driven. 

We can use this fact. For example, if we drive only bar- 
riers Wi, then state probabilities remain constant, which is 
trivial to check, and hence currents are zero on average. The 
dual counterpart of this trivial result is the no-pumping theo- 
rem lfl8l [l9ll stating that periodic driving of node energies Ei 
cannot induce a current on average. By analogy, breaking the 
detailed balance by a time-independent non-conservin g fo rce 
one can arrive at the generalized no-pumping theorem [20]. 



B. Fluctuation Relations for Currents 

When the detailed balance condition is broken, kinetic rates 
on a cyclic graph can still be parametrized by a set of N en- 
ergy parameters Ei, N barrier sizes Wi, and, in addition, by 
a set of free energy inputs on each clockwise (anti-clockwise) 
transition, /+ (f^~), such that 



The entropy production, A, sometimes is called the macro- 
scopic thermodynamic affinity [1] and cycle force £2]. This 
reversed entropy production is consistent with the observation 
that the dual and original model have the inverted energy land- 
scape. 

Recently, Sinitsyn et al have shown ll2lll that for arbitrary 
closed networks with two driving protocols, such that either 

(I) only barrier sizes are driven, or 

(II) only node's energies are driven, 

while keeping all other parameters, including A, constant 
(which is equivalent in our case to keeping either ln(fc z + / k~ +1 ) 

or hi{k d+ /k^T^ = hi(k^/kf) constant), the following ex- 
pression relates probabilities of conserved current, J c , for for- 
ward and backward driving protocols: 



Pf(J c ) 

M-J c ) 



(28) 



This Fluctuation Relation for Currents is a generalization of 
previously known relation for steady state currents [I22I to 
two classes of explicitly time dependent driving protocols. Its 
proof for graphs of arbitrary geometry in [21 1 was based on 
symmetries of the twisted master operator and its supersym- 
metric counterpart. Duality now allows us to achieve the same 
result for the cyclic graph with much less efforts. Indeed, the 
case (I) directly follows from the well known work relation 
[ElEi: 



Pf(>V) 
Pb(-W) 



(29) 



k- 



ke KE t -w i+ fn. (25) 



The work W in our case is defined during a periodic driving 
cycle, which, in the large time limit n — > 00 can be written as 
M 

W = AJ C + f dtiY^Ei (t)r)i(t)], (30) 
Jo 4 

where r/(t) is the random process given by rji(t) = 1 if the 
node i has the particle inside it at time t and otherwise zero. 

When only barriers are driven, the second term in Eq. (l30i 
is identically zero, so the work per period is given by W = 
AJ C , just like in the static case. Since work in case (I) is 
proportional to J c with a constant factor, A, the relation ( l28| > 
for J c is consequently satisfied. 

The case (II) does not straightforwardly follow from 
Eq. (f30b , however, using the previously discussed duality be- 
tween protocols with only energies and only barriers driven, 
we should conclude that Eq. ( |28T i should be satisfied for pro- 
tocols with only energies driven, which concludes its proof. 
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VI. DUALITY IN TWISTED FOKKER-PLANCK 
EQUATION 



d x )e-^(- X + d x )e*, and L* = ( X ~ $ - d x )e*( X 
tp' — d x )e~^ , which lead to 



Here we will explore the duality relations that are hidden 
in the twisted Fokker-Planck equations. The additional goal 
is to obtain the continuous limit of the recently discovered 
fluctuation relations for currents I12lll which so far have been 
discussed only in the framework of the discrete graph models. 

Let us take the jump (transition) spacing Ax, denote 
x/Ax — »■ x an d consider the continuous limit: Ax — > 0. 
The twisted Fokker-Plank equation is thus defined as the con- 
tinuous limit of the discrete twisted master equation, which 
reads 

d t p x (x,t) — L x p x (x,t), with (31) 
L x = { X -d x )[A{x,t) + {x-dx)B{x,t% 

where L x is the twisted Fokker-Planck operator with auxil- 
iary counting parameter \- Px( x ^) ' s tne generating func- 
tion of probabilities of a particle being at position x at time 
t having produced a current J by this time. A(x, t) — 
Ax[k + (x, t) — k~(x, t)] denotes the drift term and B(x, t) — 
Ax 2 [k + (x 1 t) + k~(x,t)]/2 indicates the diffusion. This 
twisted Fokker-Planck equation, as the continuous version 
of the twisted master equation for cyclic graphs, describes 
the one-dimensional stochastic motion with periodic bound- 
ary condition. 

Let us define two auxiliary functions ll25l l26tl : ip(x) — 
— f£ dyA(y)/B(y), and <p(x) — lnB(x) + ip(x), such that 
the above equation is characterized equally well by these two 
functions as by the drift, A, and the diffusion, B. These two 
auxiliary functions render us able to recast Eq. (13D as 

d tPx (x,t) = ( X -d x )e-^(x-d x )efp x (x,t). (32) 

Further, by introducing two twisted operators: Q x = (x — 9 x ) 
and J x = e~^'( X — d x )e v , we can reconstruct L x , as well as 
the other three operators as follows, 



QxJx 



( X -d x )e-^(x-d x )e^; 



(33) 



fd _at;t 
^-x ~ ^x J x 



L -x - J x Q-x 



(- x -d x )e*{- X -d x )e^- (34) 
r ' VC = e^(~x-d x )e^(~x~d x ); (35) 



L x = J x Q x = e-V(x-d x )e'<>(x-d x ). (36) 

Following the same arguments as for the discrete master equa- 
tion, it is straightforward to find that these four operators have 
the identical eigenvalues as well and share the relations de- 
picted in Fig. |2] A continuous analog of Eq. ( fTSI l thus is given 

by 



^(x, t) = t), ip d (x, t) = -<p(x, t), 



(37) 



The relations $1% and d20l > also hold true in the contin- 
uum limit, if we define the evolution operators for the twisted 
Fokker-Planck operators accordingly. We further note that 
Eq. (l35l l and d36l l are equivalent to L_ x = (~x + f' ~ 



L-x — L-x+<p'i anc ^ L x — L x _^,, (38) 
with the prime meaning the derivative with respect to x. 



A. Fluctuation Relations for Currents in one dimension 

Duality and symmetry relations on the level of twisted op- 
erators allows us to derive various relations between statistics 
of currents in original and dual models driven by periodic, 
and otherwise arbitrary, protocols. Considering the isospectral 
property of these twisted operators, their cumulant generating 
functions satisfy: 



(39) 



with the dual relations: 



i/4(t) = -Mt) = -Mr - 1) = - 1), 

<4 (i) = - # (t) = -Vb (t - t) = 4 (r - 1). (40) 

After the inverse Legendre transformation lUol IT5I and un- 
der some constraint, we are capable to obtain various fluctu- 
ation relations for currents. For instance, the second equal- 
ity in Eq. (|39l leads to the identical distribution of opposite 
currents generated by original and dual models with time- 
reversed driving protocols: 



Pf(J c ) = p£(-J c ), 



(41) 



which is the same as the Eq. (|23l obtained for discrete graph. 
We note that Eq. ( |4TT > holds true for arbitrary time-dependent 
protocols. Now, let us have either time-independent ip'(x) or 
time-independent ip'(x), which is similar to the two required 
protocols for obtaining Eq. (|28l in the discrete version. In this 
way, the first or the last equality in Eq. ((39) provides the rela- 
tions for current distributions generated by forward and back- 
ward protocols in the original and dual models, respectively: 



Pb(~J c ) 



„AJ C 



pj(J c ) 
pi(-J c 



-AJ a 



(42) 



Here, the constant, A, which is defined as 

§ dxA(x,t) I B(x,t), is also given by 

A — f+SW—fttfM, (43) 

which is similar to its discrete version (127V In fact, as 

\k+(x,t) - k~(x,t)\ <C k+(x,t),k-(x,t), 

A(x,t) 2[k+(x,t)-k-(x,t)] fc+(M) 

B(x,t) X k+{x,t) + k~(x,t) ^ n /c-(x,i)' 

so that, after a closed cycle, 



A(x,t) 
B(x,t) 



dx In 



Kit) 



(45) 
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It is clear that only the nonconservative part of 
A(x,t)/B(x,t) contributes to A. It plays the role of 
nonequilibrium source, which is related to the macroscopic 
thermodynamic affinity J3] and cycle force |0]. If A is zero, 
the system is at the situation that no directed transport occurs 
at the steady state lf27ll . 

B. No-Go Theorems for One-Dimensional Brownian Motor 

The duality relations, which are hidden in the twisted 
Fokker-Planck equations, shed a new light on some of 
the previously known no-pumping (no-go) theorems for a 
one-dimension Brownian motor Il26ll . If f{x) holds time- 
independent, the system has the steady state solution: p{x) ~ 
e -<p(x) so that the average current vanishes, no matter 
how arbitrarily we drive ip(x,t). Therefore, the dual re- 
lation between ip and rp offers us the no-pumping theo- 
rem stating that a current on average can not be pumped, 
also in case when ip(x) is fixed, which is the same as to 
saying that there is no-pumping situation when the ratio 
A(x, t)/B(x,t) is time-independent ll26ll . One concrete ex- 
ample is the overdamped Brownian motor with non-constant 
friction coefficient: d t p{x,t) = d x (U'(x)/-f(x,t)p(x,t)) + 
d^(kBT/j(x, t)p(x, t)). No matter how we vary the friction 
coefficient j(x, t), the constant ip(x) = — J? dyU'(y)/kBT 
guarantees the result that the average current is zero, which is 
consistent with the previously found no-go theorem ll28ll . 

Actually, under the detailed balance condition, we can write 

k+ ^ e E(x,t)-W(x+Ax,t) an( j ^ e E(x,t)-W{x-Ax,t) p Q i_ 

lowing these expressions, it is straightforward to show that 

A W'e E - w and B ~ e E ~ w , so that ip ~ W, ip ~ E. 

Thus, the duality between ip and tp is similar to the duality 
between energy wells and barriers in discrete graphs. Consid- 
ering the energy landscape between original and dual model 
is reverted, the signs before A in fluctuation relations Eq. (l42l 
are opposite. 

If after inverting the energy landscape or exchanging the 
energy wells and barriers, a model is invariant, we say this 
model is self-dual. In this case, the relation (|4TT > changes to 



Pf(J c ) = Pb(— J c ). Once additionally, the driving proto- 
col follows the time inversion invariance, i.e., B = F, we 
have Pp(J c ) = Pp(—J c ), which leads to the no-pumping 
(J c ) = naturally. In fact, this self-duality, combined 
with the time inversion invariance, recovers the supersymme- 
try criterion for vanishing currents in Brownian motors [l28ll : 
-V(x, / (t)) = V(x + L/2, f(-t)), where the driving has 
time-reversal symmetry: f{t) = f(—i) and the potential is 
self-dual: ~V(x) = V(x + L/2) with L/2 determined by its 
periodicity L. Our relation Pp(J c ) = Pp(—J c ) offers more 
information for this case because it implies that not only the 
average current but also all the odd cumulants of the current 
are zero. 

VII. SUMMARY 

Dualities have proven to be a fruitful guide to many funda- 
mental problems in physics. Our duality relations for stochas- 
tic currents show that fluctuation theorems for currents exist 
beyond standard fluctuation theorems for work and entropy 
production. Existence of duality relation also provides new 
means to search for new exact results in mesoscopic statistical 
mechanics because some of already known exact results may 
have their dual counterparts, as for the case of no-pumping 
theorems. We hope our results will shed new light on the 
study of time dependent nonequilibrium transport, provide 
new insights to understanding working mechanisms of bio- 
logical molecular motors, and will be applied to optimal de- 
sign of various artificial molecular motors I129ll and nanoscale 
pumps. 
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